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Many studies have analyzed how variabihty in reproductive success affects fitness. However, each 
study tends to focus on a particular problem, leaving unclear the overall structure of variability 
in populations. This fractured conceptual framework often causes particular applications to be in- 
complete or improperly analyzed. In this paper, I present a concise introduction to the two key 
aspects of the theory. First, all measures of fitness ultimately arise from the relative comparison of 
the reproductive success of individuals or genotypes with the average reproductive success in the 
population. That relative measure creates a diminishing relation between reproductive success and 
fitness. Diminishing returns reduce fitness in proportion to variability in reproductive success. The 
relative measurement of success also induces a frequency dependence that favors rare types. Second, 
variability in populations has a hierarchical structure. Variable success in different traits of an indi- 
vidual affects that individual's variation in reproduction. Correlation between different individuals' 
reproduction affects variation in the aggregate success of particular alleles across the population. 
One must consider the hierarchical structure of variability in relation to different consequences of 
temporal, spatial, and developmental variability. Although a complete analysis of variability has 
many separate parts, this simple framework allows one to see the structure of the whole and to place 
particular problems in their proper relation to the general theory. The biological understanding of 
relative success and the hierarchical structure of variability in populations may also contribute to a 
deeper economic theory of returns under uncertainty}^ 



I. INTRODUCTION 

Natural selection favors traits that enhance fitness. 
But how does one measure fitness? Several studies have 
shown that it is not just average reproductive success 
that matters. Variation in reproductive success also plays 
an important role in determining long-term evolutionary 
trends. To understand the basic notions of fitness and 
evolutionary change by natural selection, one must un- 
derstand the particular consequences of different kinds of 
variation. 

The literature on variation splits into two groups. On 
one side, bits of folk wisdom dominate thinking. The 
slogan that natural selection maximizes geometric mean 
fitness is one example. Such folk wisdom is true in special 
cases. But as a guiding principle, the simple geometric 
mean slogan misleads as often as it helps. 

On the other side, a technically demanding specialist 
literature divides into numerous distinct ways of fram- 
ing the problem. Each technical expression emphasizes a 
particular aspect of variation, refining unique examples 
at the expense of providing a coherent view of the whole. 

This paper provides a tutorial on the different kinds of 
variation and their evolutionary consequences. I empha- 
size simple examples to develop understanding of tem- 
poral, spatial, developmental, and trait variation. Each 
type of variation was originally studied as a separate 
problem. In this tutorial, I follow Frank and Slatkin [Tj, 
who showed that these seemingly different types of vari- 
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ation can be understood in a unified way. The unified 
framework arises from two steps. 

First, it is relative reproduction that matters. Only 
those traits associated with relatively greater success 
than average increase over time. Relative measures of 
success induce diminishing returns: a doubling of repro- 
duction provides less than a doubling of relative success 
PQ[2]. With diminishing returns, increasing variation in 
reproductive success reduces fitness. 

Second, the different types of variation can be ex- 
pressed as different levels within a unified hierarchy [T]. 
Variable success in different traits of an individual affects 
that individual's variation in reproduction. Correlation 
between different individuals' reproduction affects varia- 
tion in the aggregate success of particular alleles across 
the population. Temporal, spatial, and developmental 
variation affect the way in which individual variations 
combine to determine the overall variability in the num- 
ber of copies produced by a particular allele. 

I also discuss the relation of economic theories of risk 
and uncertainty to evolutionary theories of variability. 



II. RELATIVE SUCCESS INDUCES 
DIMINISHING RETURNS 

The success of genes and of traits must ultimately be 
measured by their relative frequency in a population. 
The calculation of relative frequency leads to surprising 
consequences when there is variability ^{5]. 

To illustrate the problem, consider two alternative 
types in a population, Ai and A2 . I will often refer to the 
alternative types as alleles at a genetic locus. However, 
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Box 1. Topics in the theory of natural selection 

This article is the first in a series on natural selection. 
Although the theory of natural selection is simple, it remains 
endlessly contentious and difficult to apply. My goal is to 
make more accessible the concepts that are so important, yet 
either mostly unknown or widely misunderstood. I write in 
a nontechnical style, showing the key equations and results 
rather than providing full derivations or discussions of math- 
ematical problems. Boxes list technical issues and brief sum- 
maries of the literature. 



the same analysis would apply to any heritable alterna- 
tive types in a population that have the same essential 
properties as alleles. 

Some simple notation helps to express the argument. 
Each definition uses subscripts to associate with the al- 
ternative alleles, Ai and A2, respectively. Define qi and 
q2 as the allele frequencies, such that qi + q2 — 1. Let 
Ri and R2 measure reproductive success, the average 
number of descendant copies produced by each parental 
allele. The average reproductive success in the popu- 
lation is R = qiRi + (72^2- The success of individual 
parental copies has a random component. Thus, all of 
the measures of reproductive success fluctuate randomly. 
Throughout, the unqualified words average and mean re- 
fer to the arithmetic average. 

The frequency of Ai after one round of reproduction 



q[=qi{Ri/R)^qiFi 



(1) 



where F measures relative success. I use the word fitness 
for relative success. This equation shows that fitness ul- 
timately determines the success of an allele. Reproduc- 
tive success, R, infiuences fitness. But the key relation- 
ship between reproductive success and fitness is mediated 
through the definition for fitness 



Fi = Ri/R 



Ri 



qiRi + 92^2 



(2) 



Fig. [T] illustrates the two key properties of fitness. 
First, fitness increases at a diminishing rate with a rise in 
reproductive success [HE]- Put another way, the fact that 
fitness is a relative measure means that linear changes in 
reproductive success translate into nonlinear changes in 
fitness. 

Second, the curvature of the relation between repro- 
ductive success and fitness is frequency dependent [T] . A 
rare type has a nearly linear relation between reproduc- 
tion and fitness. A common type has a very strongly 
diminishing relation between reproduction and fitness. 
This means that rare and common types are influenced 
differently by the consequences of variability, because 
more strongly diminishing returns cause variability in re- 
productive success to impose a greater penalty on fitness. 

Fig. [2] shows that diminishing returns cause a loss of 
fitness. In the figure, expected reproductive success, R, is 
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FIG. 1. The curvature of fitness versus reproductive success 
depends on allele frequency. The plots here illustrate Eq. ^ . 
The numbers above each curve show gi , the frequency of the 
allele for which the relationship is plotted. If gi = 0.1, the 
relationship follows the upper curve for allele Ai; the lower 
curve can then be interpreted as the relationship for allele A2 
with frequency 52 = 0.9. The difference between the upper 
and lower curves illustrates the frequency dependence of the 
relation between fitness and reproductive success. Note that 
there is little curvature when an allele is rare, which leads to 
an advantage for rare types. Redrawn from Frank and Slatkin 
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Reproductive success, i?j 

FIG. 2. Increasing variation in reproductive success reduces 
fitness. Expected reproductive success is fj,. Fluctuations of 
±(5 occur. Positive fluctuations return a smaller gain in flt- 
ness than the loss sufi'ered from a negative fluctuation. Thus, 
equally frequent positive and negative fluctuations return a 
net loss. Redrawn from Frank and Slatkin J^. 



/i. Deviations of ±(5 occur, with increases and decreases 
at equal frequencies. The gain in relative fitness, F , for 
an increase of 5 units of reproductive success is less than 
the corresponding loss in fitness when reproductive suc- 
cess is reduced by 5. Expected fitness therefore declines 
as the frequency and magnitude of deviations increase. 
Note that the discount to fitness rises as the curvature 
between fitness and reproductive success increases. 
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III. REPRODUCTION MULTIPLIES AND 
VARIATION REDUCES SUCCESS 

Suppose an individual has two offspring. Each of those 
offspring has two offspring. The original individual has 
four grandchildren. Compare that output to a second 
individual that has three offspring, and each of those off- 
spring has one offspring. The second individual has three 
grandchildren. In each case, the average reproduction per 
generation is two offspring. However, the individual with 
less variable reproduction has greater success than the 
individual with more variable reproduction. 

The difference occurs because reproduction multiplies 
over time. The value 2 x 2 = 4 is greater than 3x1 = 
3, even though the arithmetic averages are the same in 
each case. In multiplicative series, variation reduces the 
multiplicative product. Rather than measuring success 
by the arithmetic average, such as (3 + l)/2 = 2, the 
proper average in a multiplicative series is the number 
that, when used to multiply in each generation, gives 
the total output. In the reproductive series over two 
generations of 3 x 1 = 3, we need a number that when 
multiplied by itself gives three. This works out as -s/S x 
•\/3 = 3, so the multiplicative mean is 73 « 1.73. The 
multiplicative mean is usually called the geometric mean. 

A. Approximation for the geometric mean 

A simple approximation of the geometric mean is often 
useful. Suppose the arithmetic mean of a series is /i, and 
the variance of the series is . Then the geometric mean 
is approximately fi — cr^/2/i. This approximation shows 
clearly how the variance reduces the geometric mean. 

For example, in the series 3x1 = 3, the mean is 2, 
the variance is 1, the true geometric mean is 1.73, and 
the approximation gives the geometric mean as 1.75. The 
smaller the variance is in relation to the mean, the better 
the approximation. In evolutionary models, one often 
makes the assumption that average reproductive success 
is close to one, ~ 1, so that the approximation for the 
geometric mean often appears as ^ — (t^/2. 

B. The geometric mean principle 

Reproduction multiplies. The greatest multiplica- 
tive series has the highest geometric mean reproduction. 
Thus, the type with greatest long term fitness would ap- 
pear to be the type with the highest geometric mean 
reproductive success. That conclusion is often called the 
geometric mean principle. 

The highest geometric mean is sometimes associated 
with greatest evolutionary success. However, geometric 
mean reproductive success by itself often misleads, be- 
cause it hides more than it helps. The problem of evo- 
lutionary success and fitness turns out to be more subtle 
and more interesting. The following sections explain why. 



IV. POPULATION REGULATION AND 
RELATIVE SUCCESS 

The total resources available to the population limit 
reproductive success. That density dependent competi- 
tion causes the reproductive success of each type to be 
influenced by the reproduction of other types. For that 
reason, one cannot simply multiply the reproductive suc- 
cesses of each type independently and then compare the 
long-term geometric means. Instead, each bout of den- 
sity dependent competition causes interactions between 
the competing types. Those interactions depend on fre- 
quency. Reproduction of a rare type has little compet- 
itive effect on a common type. Reproduction of a com- 
mon type has a strong competitive effect on a rare type 
(Fig.[l])- 

The fitness measure of relative success in Eq. ^ ac- 
counts for density dependent interactions. The particular 
way in which density dependent competition arises has 
important consequences. 



V. EXPECTED CHANGE IN FREQUENCY 

We can circumscribe the main conceptual issues by 
focusing on the expected (average) change in allele fre- 
quency. Each process can be studied with respect to 
whether it tends to increase or decrease the expected fre- 
quency change. 

Eq. ([T]) gives the definition for the allele frequency in 
the following generation, q'l — qiFi. The change in allele 
frequency is Agi — q'l — qi = qi{Fi — 1). Using the 
definition for fitness in Eq. ([2]) allows us to write the 
change in frequency as 

A91 = ,^<l. (^^) ■ (3) 

The reproductive successes fiuctuate randomly. Because 
those random fluctuations occur in both the numerator 
and the denominator, there is no simple way to express 
the exact change in frequency. If we assume that the 
fluctuations in success are small relative to the average 
success, and we normalize all of the successes so that 
they are close to one, then we can write the approximate 
expected change in frequency as 

E(A(?i) w gi92 {(mi - M2) + [cov(i?2,i?) - cov(i?i,i?)] } , 

(4) 

where and /i2 are the expected reproductive successes 
for types 1 and 2, respectively. This equation, from Frank 
and Slatkin [Ij, is equivalent to an approximation given 
by Gillespie [5]. 

The expected change is the average tendency. Because 
of the inherent fluctuations in success, the actual change 
in frequency in any generation may be in the opposite 
direction of the expected change. Over long time periods, 
three different patterns of evolutionary dynamics occur. 
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First, if the random fluctuations in the average repro- 
ductive successes of the types are large relative to the 
directional bias in Eq. Q, then randomness dominates. 
The type frequencies will bounce up and down in a nearly 
neutral way. Eventually, one type will become fixed and 
the other will disappear from the population. If we start 
at frequency gi, then the probability that type Ai fixes is 
qi , and the probability that type A2 fixes is (72 = 1 — . 
True fixation only occurs in finite populations. In infinite 
populations, the related notion of quasi-fixation arises as 
the frequency of a type becomes very close to one. To 
keep things simple, I ignore important technical distinc- 
tions between finite and infinite populations ^ [5] . 

Second, if the directional bias is much larger than the 
fiuctuations in the average reproductive successes of the 
types, then the type frequencies change in an almost de- 
terministic way. If the direction of change remains the 
same across the range of type frequencies, then the type 
with the greater expected success will usually become 
fixed in a relatively short period of time. If frequency 
dependence provides a sufficient advantage to the rare 
type, then the direction of change may shift with type 
frequency in a way that tends to maintain both types in 
the population. 

Third, if random fluctuations are of roughly the 
same magnitude as the directional bias, then frequency 
changes combine both directional and random aspects. 
In some cases, frequencies will fluctuate, and both types 
will remain in the population for a very long time. In 
other cases, frequencies may fluctuate for a significant 
period of time, but eventually one type or the other will 
become fixed. Fixation will be biased toward the type 
favored by the directional tendency set by Eq. Q . How- 
ever, the other type may occasionally fix because of the 
random fluctuations. 

I emphasize the major processes that influence the di- 
rectional tendency in Eq. Q. In particular, the hier- 
archical structure of variability sets the directional ten- 
dency, which in turn shapes the qualitative patterns of 
dynamics. I make only brief comments on long-term evo- 
lutionary dynamics, which require technical analysis of 
mathematical models to evaluate fully [H [S] . As noted 
in the previous paragraphs, long-term dynamics include 
issues such as the probability that a particular genotype 
becomes fixed and the maintenance of polymorphism. 



For example, the reproductive success i?i is the average 
success taken over all copies of the allele Ai. Similarly, 
i?2 is the average over all copies of ^2- To analyze the 
variability in the aggregate success of allele Ai, one must 
consider the variability in the success of each copy of Ai 
and the correlations in success between different copies. 
The analysis for A2 must consider the variability in the 
success of each allelic copy and the correlations in success 
between the alleles. 

We can relate the variances of individual reproductive 
success to the variances and covariances for the different 
allelic types by writing 

var(iJi) = piaf (5a) 

var(i?2) = /02O-2 (5b) 

C0v(i?i,_R2) = Pl2CTl(T2 (5c) 

where pi, p2, and pi2 are the correlations in reproductive 
success between randomly chosen pairs of Ai, A2, or Ai 
and A2 individuals, respectively. 

Frank and Slatkin jlj introduced this explicit parti- 
tioning of variances and covariances for types into their 
individual components. Any realistic analysis of variabil- 
ity must make explicit the individual level fluctuations 
and the associations between individuals. Although this 
explicit treatment of variability is fundamental, the par- 
titioning of variances and covariances in Eq. ([s]) has often 
been regarded as some sort of highly technical or special- 
ized analysis. This mistake has limited progress in un- 
derstanding fitness with respect to spatial and temporal 
fluctuations in success. 

If we combine Eq. Q and Eq. ([5|, and keep things 
simple by assuming that the correlation between types is 
zero, pi2 = 0, we obtain 

E(A(7i) w qiq2 {{pi - qiPicrl) - {p2 - q2P2C^D} , (6) 

which means that, on average, type Ai increases in fre- 
quency when 

Pi - qiPicrj > P2- q2P2crl- (7) 

The following sections show that different kinds of vari- 
ability can be understood by the hierarchical partition- 
ing of associations between traits within an individual 
and associations between different individuals. Variabil- 
ity interacts with the processes of density dependent pop- 
ulation regulation. 



VI. HIERARCHICAL STRUCTURE OF 
VARIABILITY 

Populations have a naturally nested hierarchical struc- 
ture when considering genetics. Populations are com- 
posed of genotypes, and genotypes are composed of indi- 
viduals. In a haploid model with two alternative alleles, 
there are two genotypes in the population, and numerous 
copies (individuals) of each allelic type. The hierarchical 
structure of variability makes explicit the variances and 
the correlations at different levels in the hierarchy. 



VII. TEMPORAL VARIABILITY 

Dempster [6] introduced a model of temporal variation 
in which all alleles of the same type have the identical 
reproductive success within a generation, pi = p2 = I, 
and there is no correlation between types, pi2 — 0. In 
this haploid model, each individual has one allele, either 
Ai or A2 ■ The condition for the expected increase of type 
1 from Eq. Q is 

Pi - qi<yl > P2- q2<yl- (8) 
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This equation illustrates the rare-type advantage induced 
by density dependent population regulation. When the 
frequencies are equal, 9i = 92 = 1/2, then the condition 
favors the type with the higher geometric mean fitness, 
/ii — cri/2 versus ^2 — ct|/2. 

As the frequencies approach one boundary, gi — > 
and (72 ^ 1, the condition to favor Ai becomes /^i > 
/i2 — cr|- At the other boundary, the condition favoring 
Ai becomes ^i—af > fi2- Thus, the directional tendency 
often shifts with frequency. 

In spite of the inherent rare type advantage, polymor- 
phism is not maintained in this haploid model [THS]. The 
high variance in fluctuations eventually causes one of the 
types to fix (or to become nearly fixed in an infinite popu- 
lation). The type with the higher geometric mean success 
has the advantage at the frequency midpoint. That type 
fixes with higher probability. If the geometric means for 
the two types are close to each other, then frequencies 
may fluctuate for a long time, and the bias toward fix- 
ing the favored type is relatively weak. If the geometric 
means for the two types are significantly different, then 
fixation happens sooner and with a stronger bias toward 
the favored type. 



VIII. CORRELATIONS AND GENOTYPIC 
HOMEOSTASIS 

In the previous section, I assumed that all individu- 
als carrying the same allele have the same reproductive 
success in each generation. In that case, all of the varia- 
tion arises from the response of an allele to environmental 
fluctuations, with no variation between individuals of the 
same genotype. No variation means that individuals of 
the same type are perfectly correlated, pi — p2 = 

Alternatively, different individuals of the same type 
may respond differently to environmental fluctuations. 
There are many ways to express individual variation. For 
example, individual responses may fluctuate about the 
long term arithmetic mean, /z, and the pairwise corre- 
lation between individuals in each generation may be p 
[1]. In that case, the variance in the average reproductive 
success of Ai is pitr^, with a similar expression for ^2- 

Reduced correlation between individuals lowers the 
variation in the average success of a type. That rela- 
tion arises from the fact that the variance of an average 
is reduced by the number of uncorrelated observations 
in the sample. We can express the effective sample size 
of uncorrelated observations as n* — 1/p, so that the 
variance of the mean, a'^/n* , is pa^. 

One can think of the pairwise correlations between in- 
dividuals of the same genotype as the genotypic home- 
ostasis. If all individuals of a genotype respond in ex- 
actly the same way to each environmental state, then 
the correlation between pairs of individuals is perfect, 
p = 1. That perfect correlation increases the variance 
in the average reproductive success of the genotype. One 
can think of such strong correlation as strong homeostasis 



or canalization of development for the genotype. By con- 
trast, weak correlation between individuals of the same 
genotype, with low values of p, corresponds to greater de- 
velopmental fluctuations and relatively weaker genotypic 
homeostasis. 

Given the variances in the average reproductive suc- 
cesses of the types as pcr^, the condition for Ai to be 
favored is given in Eq. ([7|. In a haploid model, the long- 
term bias in fixation depends on the relative geometric 
means derived when frequencies are equal, qi — q2 = 1/2, 
yielding the condition for Ai to be favored as 

fii - pial/2 > P2- P2(yl/'2- (9) 

This expression shows that temporal fluctuations favor 
reduced genotypic homeostasis, with low values of p. A 
type with low average reproductive success, p, can be 
favored if it also has low genotypic homeostasis, p, re- 
ducing its variance in average reproductive success suffi- 
ciently to give it a higher geometric mean fitness than its 
competitor. Reduced genotypic homeostasis is a general 
expression of the widely discussed problem of bet hedging 
(see Box [2]). 



IX. DEVELOPMENTAL VARIABILITY 

Gillespie [TO] introduced a model in which the repro- 
ductive success of each of the N haploid individuals in 
the population depends on its interactions with the en- 
vironment during development. The reproductive suc- 
cesses of different individuals are independent because, 
by Gillespie's assumptions, difl^erent individuals experi- 
ence different conditions and develop in an uncorrelated 
way. Nevertheless, the finite population size ensures that 
an individual's reproductive success correlates with the 
average reproductive success of its genotype. 

The correlation of two randomly chosen Ai alleles is 
Pi — \/ {Nqi ) , because there are Nqi individuals of type 
Ai, and hence a chance l/{Nqi) of choosing the same 
individual twice. By the same reasoning, p2 = l/{Nq2). 
The correlation between types is zero, because different 
individuals experience different conditions. Substituting 
these values into Eq. ([s]), wc find that Ai increases for 
any allele frequency when 

pi - al/N >p2~ (jI/N. (10) 

Since this condition no longer depends on allele frequen- 
cies, it is sufficient to describe long-term evolutionary 
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Box 2. Optimal phenotypes in response to environ- 
mental variability 



Individuals may be able to match their phenotype to par- 
ticular environments. Phenotypic plasticity occurs when an 
organism can sense the particular environmental state and 
adjust its traits accordingly. If organisms do not adjust their 
phenotypes in response to the particular environmental state, 
then they may produce a stochastic response tuned to the 
pattern of fluctuation. A stochastic response is sometimes 
called bet hedging. 

Bet hedging can increase the aggregate success of a geno- 
type or strategy. Suppose, for example, that the environment 
is equally likely to be in one of two states. For each state, 
there is a different optimal phenotypic response. However, the 
organism cannot adjust its phenotype in response to the par- 
ticular state. If each individual of the genotype has a random 
component to its phenotypic expression, then in each gener- 
ation some individuals will match the environment with the 
best phenotype and some will not. The mixture of phenotypic 
expressions reduces the variance of the aggregate success of 
the genotype by reducing the correlation between individuals 
of that genotype. 

The concept of reduced correlation between individuals of 
a genotype highlights an essential aspect of the bet hedging 
problem. Frank and Slatkin [l] emphasized the general point 
about correlations between members of a genotype, as shown 
in Eq. ([5|. McNamara [11] independently described a similar 
interpretation, but referred to the process as kin selection, a 
label that I would avoid in this case. Correlations between 
relatives do matter, but not in the way that one usually as- 
sociates with the costs and benefits of social behavior in kin 
selection. 

Bet hedging can also arise within an individual. For ex- 
ample, the individual's alternative traits may be expressed 
stochastically in separate bouts of resource acquisition, in a 
way that reduces the overall variance in success. Reduced 
variance in resource acquisition typically provides increased 
expected reproductive success, because success rises in a di- 
minishing way with resources. 

To get started on the literature, here are a few recent 
overviews for phenotypic plasticity [12fll9j and bet hedging 
[20 '-"26 . An interesting information theory approach may pro- 
vide a useful connection between these subjects |27] . 

All of these cases analyze phenotypic adjustment or phe- 
notypic stochasticity. In these cases, one must also account 
for the diminishing relation between reproductive success and 
fitness (Fig. [T| and the ways in which the correlations between 
individuals determine the mean and variance of aggregate suc- 
cess for each type (Eq. |5|. 



advantage without the need to consider frequency de- 
pendence. An allele with a long-term advantage is more 
likely to become fixed than a neutral allele with the same 
initial gene frequency. 

Gillespie [TD] presented this model of individual devel- 
opmental variation as a separate problem from the gen- 
eral analysis of fluctuations in reproductive success. The 
analysis here, from Frank and Slatkin 1 , shows that indi- 
vidual variation is just a special case of the general model 



of temporal variation. One obtains the case of individ- 
ual variation by properly calculating the correlations in 
reproductive success between individuals. 



X. SPATIAL VARIABILITY AND LOCAL 
POPULATION REGULATION 

The classical Dempster 6] model of temporal variation 
assumes that density dependent regulation occurs in one 
large population. In that model, density regulation in- 
duces frequency dependence that favors the rare geno- 
type. I mentioned earlier that, in spite of the rare type 
advantage, one of the types eventually becomes fixed, 
because the random fluctuations in frequency are too 
strong relative to the directional tendency of evolution- 
ary change. Fixation is biased toward the type with the 
highest geometric mean. 

In a different model, Levene [55] showed that spatial 
variation does maintain genetic polymorphism. In the 
Levene model, there are many independent spatial lo- 
cations. Each location has its own independent density 
dependent competition for resources. 

Gillespie [211 HHj showed that one can think of the 
Levene model of spatial variation as the sum of K in- 
dependent models of temporal variation. If there is only 
one patch, K = 1, then reproductive successes fluctuate 
over time in that patch, and all competition occurs in 
that single patch. This model is identical to the classical 
Dempster model for temporal variation. 

As K increases, each independent patch fluctuates 
with the same rare type advantage of the classical Demp- 
ster model. The total fluctuation in each generation is 
the average of the fluctuations over all patches. Because 
the patches fluctuate independently, the variance of the 
average fluctuation over the entire population is reduced 
by 1/K. This reduction arises because the variance of 
the mean for K independent observations is the variance 
of each observation divided by K. 

Gillespie [291 ISD] provided the full analysis for this av- 
eraging over K patches. However, he did not provide a 
simple interpretation or a simple expression for how fluc- 
tuations lead to a particular level of polymorphism when 
the number of patches is large. 

Frank and Slatkin [T] noted that, as K becomes large, 
the population-wide fluctuations in each generation be- 
come small because of the averaging effect over the many 
patches. Thus, we can treat Eq. (|6| as an essentially de- 
terministic process. The rare-type frequency dependence 
now dominates. The equilibrium frequency of types can 
be obtained by solving E{Aqi) = 0, which yields 

<Zi ^ Ml - + P2gi ^^j^^ 

as given by Frank and Slatkin [1]. Here, each p is the cor- 
relation between copies of an allelic type measured within 
each patch. This result shows that naive comparison of 
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geometric mean success is not sufficient to understand 
evolutionary outcome. 



XI. TRAIT VARIABILITY WITHIN 
INDIVIDUALS 

The theory above takes an individual's average and 
variance in reproductive success as given parameters. 
However, the actual processes that lead to individual 
means and variances arise from the way that individ- 
uals acquire resources and produce offspring, including 
acquisition of food, protection from predators, and so 
on. To analyze the full hierarchy of causes for variabil- 
ity, we should begin with the question: How does the 
allocation of an individual's resources among alternative 
traits influence that individual's mean and variance in 
reproductive success? 

I use the example of traits for resource acquisition. 
The same analysis applies to any trait that influences 
reproductive success, such as defense against parasites. 



A. One trait 

Let us start with a single trait. The return of resources 
on investment has a random component, S. The ran- 
dom component of resource acquisition affects reproduc- 
tive success by an amount f{S). Then a simple way to 
write the reproductive success is 



(12) 



If we assume that the random fluctuations, S, have a 
mean of zero and a variance of Vx , and that the fluctua- 
tions are relatively small, then the average reproductive 
success is approximately 



M « 1 + f"Vj2, 



(13) 



where /" is the second derivative of / evaluated at zero 
[3T1[32]. Typically, one assumes that fluctuations in traits 
for resource acquisition have a diminishing return shape 
as in Fig. [2] in which case /" < 0. Thus, greater fluctua- 
tions, Vx, reduce expected reproductive success. AU else 
equal, resource acquisition strategies with less variabil- 
ity yield higher average reproductive success than those 
strategies with more variability. 

The variance in an individual's reproductive success is 
approximately 



« var(/'5) = /'V,, 



(14) 



where /' is the flrst derivative of / evaluated at zero. 

A full evaluation of fitness requires specifying the 
means, variances, and correlations between all individ- 
uals in the population. The correlations must be evalu- 
ated in relation to the heritable types we are following 
over time. The earlier sections provided the methods for 
studying evolutionary dynamics in relation to fitness. 



Here, to keep the focus on trait variability within in- 
dividuals, I give only the geometric mean reproductive 
success for an individual. The geometric mean for each 
individual accounts for the average and variance in indi- 
vidual reproductive success but neither the correlations 
between types nor the role of density dependence in fit- 
ness. Assuming that the fluctuations in returns, V^, are 
relatively small, the geometric mean reproductive success 
is approximately 



G 



aV2M « 1 + (/" - r ) VJ2. (15) 



B. Two traits 

How should an individual invest in two different traits 
that provide additive returns? Let reproductive success 
be 



R^x[l + fiS)]+y[il~^)+gie)] 



(16) 



for investment amounts x + y — 1, with x and y the frac- 
tions of total resources invested in each trait. Here, 7 is 
the discount in expected return for the second trait, and 
e is the random fluctuation associated with the second 
trait. The discount, 7, and the fluctuation, e, are small 
relative to the baseline return of one. The mean of e is 
zero, and the variance is Vy. I assume that 6 and e are 
uncorrelated. Box [3] fists some of the intermediate steps. 
Here, I focus on the key result. 

If we assume that the random component of each trait 
is the same, Vx = Vy, and f = then the geometric 
mean is 



/z-(TV2/i«G + B(x,2/), 



(17) 



where G is the geometric mean given in Eq. ( 15 ) for allo- 
cating all resources to the first trait, a; = 1, and B{x,y) is 
the benefit obtained when mixing allocation of resources 
between the two traits such that x + y = 1, with 

y) = r [1 - {x^ + J/')] Vx/2 - y-1. (18) 

If we optimize B to obtain the best mixture of allocations 
between the two traits, we obtain 



* 1 ^ 



(19a) 
(19b) 



where 7 is the dis cou nt in expected return for the second 
trait given in Eq. ( 16 1, and cr^ is the variance in individual 



reproductive success per trait given in Eq. ( 14 ) 
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Box 3. Trait variability with two traits 

This box shows the details that lead from Eq. (|16|l to 



Eq. (17 1. Starting with Eq. (16 1 for R, the average reproduc- 

(20) 



five success is approximately 

M ^ 1 - y7 + ^f"V./2 + yg"Vy/2 
The variance in success is approximately 



(T^ ~ va.r{x f'S + yg'e) 
= x'^f'^V^ + y'^g'^Vy. 

The geometric mean is approximately 

M - aV2/i ^ 1 - 2/7 + ^ (/" - V./2 + y {g" - yg'^) Vy/2. 

If we assume that the random component of the two traits is 
the same, Vx = Vy, and f = g, then the geometric mean is 

M - aV2M ^ 1 - y7 + /"K/2 - {x^ + 2/')/'V./2 
= 1 + (/" - n K/2 

+ [1 - (x' + y')] K/2 - y-y. 



The substitutions given in the text lead directly to Eq. ( 171 



It pays to invest some resources in the trait with lower 
expected return as long as < 1. The lower ex- 

pected return is offset by the reduced variance obtained 
from averaging the returns over two uncorrelated traits. 
In both biology and financial investing, returns tend to 
multiply over time. Thus, reduced fluctuations enhance 
the multiplicative (geometric) average return. In finan- 
cial investing, the central role of the geometric mean is 
well known in theory [33] . but often ignored in practice 



In this case, I have assumed 7 = 0.02 and cr^ = 0.1. 
From Eq. (19), the optimal allocation to the two traits 
is X* = 0.6 and y* — 0.4. If the individual devotes a 
fraction 0.6 of its investment to the first behavior and 
a fraction 0.4 of its investment to the second behavior, 
then it obtains an average return of 



a = 0.6^1 +0.4^2 = 0.992. 



(21) 



The variance in return is obtained by noting that, when 
one splits allocation between two uncorrelated returns, 
Ri and R2, each with variance cr^, the variance is 



b = va,r{xRi + j/i?2) = [x^ + y^]a^. 



(22) 



Using the optimal split 0.6 versus 0.4 for x and y, and 
the value cr^ = 0.1 above, the variance is 6 = 0.052. The 
geometric mean is now approximately 



a-b/2^ 0.992 - 0.052/2 = 0.966. 



(23) 



This mixture of behaviors therefore returns a higher geo- 
metric mean of 0.966 than when all investment is devoted 
to the higher yielding first behavior, which has a geomet- 
ric mean of 0.95, or when all investment is devoted to the 
lower yielding behavior, which has a geometric mean of 
0.93. This example illustrates the benefit of diversifica- 
tion when success multiplies over time. 

Note that reproductive returns are linear in this ex- 
ample. The entire benefit of diversification arises from 
the multiplicative nature of long-term success, which dis- 
counts variance. 



D. The limitation of using individual geometric 
mean success 



C. An example 

The concepts in the previous section are simple. The 
variance of an average declines with additional uncorre- 
lated components. Reduced variance provides a benefit 
when success multiplies over time. The technical expres- 
sions of those results may obscure the simplicity of the 
concepts. This section provides a numerical example. 

Suppose an organism has two different behaviors by 
which it can obtain calories. To keep the problem simple, 
assume that there is a linear relation between calories and 
reproduction, f — 1. For the first behavior, the return 
is on average /ii = 1.0 calories, with a variance in return 
of = 0.1. The second behavior has a lower average 
return of 112 = 1-0 — 0.02 = 0.98 calories, with the same 
variance of ct^ = 0.1. 

If all investment is devoted to the first behavior, then 
the geometric mean success is fii — (7^/2 = 0.95. If all 
investment is devoted to the second behavior, then the 
geometric mean success is fj,2 — cr^/2 = 0.93. 



I used the individual's geometric mean success in this 
section. That assumption is valid only when we are inter- 
ested in an absolute measure of an individual's long-term 
success in the absence of competition and relative com- 
parison with others. However, it often does not make 
sense to measure success independently of others. Evo- 
lutionary success depends on the relative contribution to 
the population by a heritable type. The earlier sections 
of this paper showed several different measures of suc- 
cess that arise from the temporal and spatial structure 
of competition and from the correlations in success be- 
tween different types of individuals. 

Suppose, for example, that the correlation between in- 
dividual copies of an allele is low, p — > 0, as in the model 
of developmental variation in which p = 1 /N and popula- 
tion size, N, is large. Then, from Eq. (|9|, we see that nat- 
ural selection favors the type with the highest arithmetic 
average, /i, independently of the individual variance in 
reproductive success, a^. In that case, it does not make 
sense to analyze the geometric mean of individual repro- 
ductive success. Instead, success depends almost entirely 
on the arithmetic mean return taken over the two traits. 



9 



If returns per trait are linear, /" = 0, then the arithmetic 
mean is 

xni+yn2 = x + y{l-j). (24) 

In this case, individuals will be favored to allocate all 
resources to the higher yielding trait, labeled as trait one 
in this example. 



XII. ECONOMIC THEORIES OF VARIABILITY 
AND RISK 

Economic theories of risk and uncertainty typically fo- 
cus on the absolute success of individuals or single agents. 
Relative success in economics concerns market share [3S] . 
However, there seems to be little economic theory about 
risk and uncertainty in relation to market share. Prob- 
lems of market share lead to many issues discussed in this 
paper. For example, relative success induces diminishing 
returns. The temporal and spatial scale of competition 
determines the proper measure of success. 

One must also consider the proper unit of analysis to 
measure success and dynamics. If one is interested in the 
absolute currency value accumulated by an individual in- 
vestor over a long period of time, then the individual's 
geometric mean success in return over successive inter- 
vals is often a good measure. If one is interested in an 
individual's purchasing power, then one must track the 
individual's currency valuation relative to the currency 
valuation among the population of individuals compet- 
ing for the same goods. 

If the individual has only a small fraction of the total 
pool of goods, then the individual's geometric mean re- 
turn provides a good measure of success. However, if one 
is tracking a corporation or agent that controls a large 
fraction of the total resource pool, then the correlation 
between individual and total success may have a signifi- 
cant impact on outcome. 

In some economic analyses, one is interested in be- 
haviors or strategies. For example, what is the relative 
success of those following a particular financial strategy 
in the investment markets? The answer depends in part 
on whether all individuals following the same strategy 
have highly correlated returns or uncorrelated returns. A 
high correlation in returns between individuals following 
a strategy increases the variance in the aggregate success 
of that strategy. Higher variance usually leads to lower 
long-term success. 

The frequency of the competing strategies also mat- 
ters. Relatively rare strategies have a low correlation 
with the population average level of success, providing a 
rare-type advantage. The consequences of the rare-type 
advantage depend on whether competition occurs glob- 
ally or locally over a series of isolated markets. These 
problems of relative success have received little attention 
in the study of economic competition. 



Box 4. Reviews and technical issues 



Gillespie [J and Ewens [5] provide excellent technical 
overviews of genetic theory for variable environments. Re- 
cent reviews |36l I37j and new theory continue to appear. 

Lande [39] developed a comprehensive approach to the 
theory of fluctuating selection. He emphasized an adaptive to- 
pography method arising from the key insight that expected 
fitness depends on average reproductive success minus the 
covariance of reproductive success with population mean suc- 
cess. That expression for fitness is the same as developed in 
Eq. Q, following from Gillespie [2] and Frank and Slatkin [T]. 

Lande [39] also summarized the complexities that arise in 
diploid genetic systems under fluctuating selection. In gen- 
eral, diploid and multilocus models require one to pay atten- 
tion to two issues [T]. First, how do the fluctuations con- 
tributed by different alleles combine within individuals to de- 
termine the average and variance in individual reproductive 
success? Second, how do multiple alleles per individual induce 
correlations in reproductive success between copies of alleles 
in different individuals? The specific results that I gave in the 
text concern haploid models. Diploid and multilocus models 
must account for these additional complexities. 

Rice introduced an exact expression for evolutionary 
change in stochastic models by expanding the scope and inter- 
pretation of the Price Equation. Tuljapurkar et al. [41] review 
the theory of variable environments in relation to demography 
and life history evolution. For economics theory, introductory 
microeconomics texts include overviews of the theory of risk 
and uncertainty. For investments, see Markowitz |42] . Okasha 
[43) provides an entry to the philosophical literature on the 
theory of uncertainty in evolution and economics. 



XIII. CONCLUSIONS 

Nearly all aspects of success include a variable compo- 
nent. To understand the consequences of that variability, 
one must study the hierarchical structure of traits within 
populations. Each individual has multiple traits. Each 
genotype or strategy has multiple individuals. Fluctua- 
tion in the success of particular traits has consequences 
that depend on the correlations between traits and the 
correlations between individuals. The aggregate variabil- 
ity of competing types affects relative success in ways 
that depend on density dependent competition, which 
causes diminishing returns and induces an intrinsic fre- 
quency dependence that tends to favor rare types over 
common types. 

The extensive biological theory of variability has dealt 
with particular aspects of the overall problem. But few 
analyses have set out the entire range of fiuctuations, how 
those fluctuations are structured in populations, and the 
particular nature of competition that shapes the conse- 
quences of fluctuations. By considering the structure of 
the entire problem, one obtains a richer understanding 
of biological fitness and its consequences for evolutionary 
dynamics. 
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Many of the biological problems of variability also arise 
in economics. The theoretical literature in economics 
made the first analyses of success when there is a vari- 
able component of returns. But the biological literature 
has advanced further in the analysis of variability, par- 
ticularly with respect to the importance of relative suc- 
cess and the hierarchical structure of competing types or 
strategies in populations. 
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